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Abstract 



An inverse scattering method based on an auxiliary inverse Sturm-Liouville problem 
recently proposed by Horvath and Apagyi [Mod. Phys. Lett. B 22, 2137 (2008)] is ex- 
amined in various aspects and developed further to (re)construct spherically symmetric 
fixed energy potentials of compact support realized in the three-dimensional Schrodinger 
equation. The method is generalized to obtain a family of inverse procedures character- 
ized by two parameters originating, respectively, from the Liouville transformation and 
the solution of the inverse Sturm-Liouville problem. Both parameters affect the bound 
states arising in the auxiliary inverse spectral problem and one of them enables to re- 
duce their number which is assessed by a simple method. Various solution techniques of 
the underlying moment problem are proposed including exact Cauchy matrix inversion 
method, usage of spurious bound state and assessment of the number of bound states. 
Examples include (re)productions of potentials from phase shifts known theoretically or 
derived from scattering experiments. 

PACS: 02.30.Zz, 02.60.Cb, 02.60.Pn, 03.65.Nk 
AMS subject classification: 34L25, 65L09, 81U40 

1 Introduction 

We consider the three-dimensional inverse scattering problem of the Schrodinger equation on 
the half-line at a fixed scattering energy [H For spherically symmetric potentials the partial 
wave expansion applies and the radial Schrodinger equations 



with the physical boundary conditions describe the scattering of two non-relativistic quantum 
mechanical objects (potential scattering). Presently k is fixed to a constant value and we 
assume that g(r) is compactly supported, i.e. g(r) = for r > a. Moreover we require that 





(1) 




(2) 



For this problem it can be shown that 

^i{r)=Cy+\l + o{l)), r^O, (3) 
ipi{r) = {Ji+i/2ikr) - tanSiYi+i/2{kr)) , r > a, (4) 

where the phase shifts 6i arise. 

In [3] it was shown that a subset of the fixed energy phase shifts, whose indices satisfy the 
Miinz condition, determines the m-function and thus the spectral function, and thereby the 
potential of an auxiliary inverse eigenvalue problem (i.e. an inverse Sturm-Liouville problem). 
Based on this proof a constructive method was suggested in [1] for the solution of the inverse 
scattering problem at fixed energy. 

In this paper we develop this method further by generalizing both the transformation of the 
fixed energy inverse scattering problem to the inverse spectral problem and the constructive 
inversion method of the inverse eigenvalue problem. In [4J only the case when one bound state 
is present in the auxiliary problem has been addressed. We examine the zero and the multiple 
bound state cases. In the latter one we find that it involves a highly nonlinear system of 
equations which is difficult to solve. One way to overcome the difficulties would be to reduce 
the number of bound states and to leave the nonlinear regime. Remarkably, this can be achieved 
in many cases by tuning the parameters found by the generalization of the original method. 
Also, an approximative argument is presented to assess the number of bound states present in 
the auxiliary problem. 

The paper is structured as follows: in the next section the transformation of the fixed 
energy inverse scattering problem is generalized. Section 3 contains a brief summary of the 
classical inverse Sturm-Liouville theory and the utilization thereof in the present problem. In 
section 4 solution methods of the inverse problem at different bound state levels are discussed. 
In section 5 an approximative method is given to determine the number of bound states and 
also the possibility to reduce their number is studied. Section 6 is devoted to illustrative 
applications, section 7 is left for a summary. 



2 Liouville transformation 



First the fixed energy problem is transformed to an inverse eigenvalue problem where the value 
of the m-function for some arguments is determined from the original fixed energy phase shifts. 
To this end we transform the radial Schrodinger equation ([T]) to the Liouville normal form (see 
e.g. [5J) by using a Liouville transformation. 



r — )■ xlr) 



(5) 



Rewriting the differential equation ([T]) in terms of the new independent variable x and dependent 
variable ipi{^) = f{^)~^fi{^{^)) yields 



<(x) 



x_ ^ J'{x) 
x"^ 



/(^)J 



q{r{x)) - _ a^f{x)_ _ f"{x) 
x2 /(x) /(x) _ 



x^ 



/(/ + ! 



r(x)^x2 



(6) 



where dot denotes differentiation with respect to r. To get the Liouville normal form of the 
Sturm-Liouville equation we need 



r{xYx'^ 



and 



X 



x^ 



const. 



fix) 
fix) 
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(7) 
(8) 
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The two conditions yield the unique solution 

x{r) = clogr + ci (9) 

and 

/(a:) = c2e^. (10) 

If we want to use the inverse spectral theory of the Sturm-Liouville equation on (0, oo) we need 
x{0) = +00 and x{a) = which in turn implies 

sgnc = —1, Ci = —clog a. (11) 

Without the loss of generality C2 = 1 is set and then the only remaining parameter c can be 
chosen arbitrarily maintaining the negative sign. 

In summary we have only one family of Liouville transformations reducing the radial 
Schrodinger equation ([1]) to the Liouville normal form, namely 

T 

x{r) = clog -, c < 0, (12) 
a 

i)i{x) = e~^ipi{ae^). (13) 

Thereby ([T]) transforms to 

- <(a;) + Q{x)^Pl{x) = + ^ M^), (14) 

with the auxiliary potential 

CL 2x 



g(x) = -e^(g(aef)-A;2). (15) 
The transformed equation can be viewed as 



S\Q{x)\y{x, A) = \y{x, A), S\Q{x)\ = + Qix) (16) 



given explicitly at 



1 f. V' 



also obtaining one of the two linearly independent solutions of the Sturm-Liouville equation as 



y\x,-\(l + l] I =Mx)- (18) 



(Later we show that this is an solution.) 

Note that in the original Horvath-Apagyi jl] method c = — 1 was taken implicitly. 

3 Inverse Sturm-Liouville problem 

3.1 Summary of the classical inverse Sturm-Liouville problem 

Spectral properties For the Sturm-Liouville equation 

- y'^{x,X) + Q{x)yaix,\)) = \yaix,\), x G [0, 00) (19) 
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with the initial conditions 

?/a(0, A) = sina 7^ 0, ?/^(0, A) = — cos a, (20) 

there exists [HI E] a monotone increasing function Pci(A), the spectral function, such that, for 
every f{x) G L^{0, oo) there exists in the L^(— oo, oo, pa{X)) norm sense 

F„(A) = l.i.m.„^oo / f{x)ya{x,X)dx (21) 
Jo 

and this is a unitary transformation, i.e. the Parseval formula 

/•oo /"OO 

/ \fix)\^dx= / |F„(A)|2cip„(A) (22) 

Jo J -oo 

holds (a theorem of Weyl). 

A property of the spectral function is the formula [HI [7] 

Pa(A) = — ^^A^/2 + p«(-oo) + ^^^ + o(l), A^oo. (23) 
TT sm a sm a 

If Q(x) G L^(0, oo) then the Sturm-Liouville operator is in the limit-point case at infinity 
and the Weyl- Titchmarsh m-function is defined uniquely by 

where y{x, A) is a solution belonging to the function space lv^(0, oo). 

The m-function is related to the spectral function through a certain kind of Stieltjes trans- 
form [6] 

sina — m(A) cosa f'^ dp ait) 

^ ^ - - cota+ / -pl2. (25) 



cosa + m(A) sina X — t 

An equivalent formulation is given for the solution of the Sturm-Liouville equation with the 
initial conditions 

/(0,A) = 1, /'(0,A) = /i<oo. (26) 

Prescribing h = —cot a one has [sina y''](0. A) = sin a and [smay'^]'{0, X) = —cos a thus 
smay^{x, X) = ya{x,X). Then Fa{X) = sin a f{x)y^{x, X)dx = smaF^{X) for f{x) G 
L^(0, oo) in the norm sense. The Parseval formula fl2^ yields 

p^(A) =p„(A)sin2a. (27) 

p'*(A) is related to the m-function by 

1 ^ rdp^ 

mix) -h J^^X-t' ^^^^ 
which formula can be inverted by the Stieltjes inversion (see e.g. XIV. §3. of [7]) 

P^Aa) - p^Ai) = -- lim / ' Im— -dX. (29) 

^ ^ ^ ^ ^ ^ TT £^0+ J)^^ m(A + ie)-h ^ ^ 

One can see that the two formulations are completely equivalent. We will use the latter one 
since it is traditionally used in the Gel'fand-Levitan construction discussed below. 

We note that in the original formalism of Horvath and Apagyi |1] /i = was taken implicitly. 
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Construction of the potential from the spectral function From the existence of the 
spectral function the Gel'fand-Levitan (GL) integral equation can be deduced [6]: 

= F{x,t) + K{x,t) + K{x,s)F{s,t)ds (0 < t < x), (30) 

Jo 

where the input symmetrical kernel is 

cos{V\x) cos(v^t)rfa(A) = -{F{x + t) + F{\x - t\)), (31) 

'OO ^ 
POD 

F{x)= cos{\^x) da {X), a(A) = p^(A) -p°'°(A). (32) 



p^{X) is defined as before while p^'^{X) is the spectral function for the zero potential with 
boundary conditions y{0) = 1, y'{0) = 0, i.e., 

(33) 

^ ^ \0, A<0. ^ ^ 

In the GL equation K{x, y) is the kernel of the transformation operator Tq q : — )■ 
realized as 

Tq,o/(x) = /(x) + / K{x,t)f{t)dt, TQ^oyofi{x,X) = yQ^h{x,X) (34) 

^0 

mapping the solutions of the Sturm-Liouville equation with Q = satisfying the boundary 
conditions 

/(O) = 1, /'(O) = 0, (35) 
onto the solutions with Q ^ satisfying 

/(O) = 1, /'(O) = h<oo. (36) 
The transformation kernel is connected to the potential Q{x) by 

Q{x) = 2^K{x,x). (37) 
3.2 m-function of the operator S[Q{x)\ 

First we show that for the above-defined potential Q{x), Q{x) G L^{0, oo) holds: 

/ \Q{x)\dx<— \q{r)\dr + -—< oo (38) 

Jo c Jo 2\c\ 

by equation ([2]), the fact that Jq r'^\q{r)\dr < a r\q{r)\dr and c < was employed. Then 
S'[(5(a;)] with Q{x) being the auxiliary potential is in the limit-point case. 
Next we prove that the functions {'?/'«(x)};=o,i,... are of the class L'^{0, oo): 

[ \'ipi{x)\'^dx = I e^~\ipi{ae~)\^dx = a\c\ f ^'^'•^J^ dr < oo, (39) 
Jo Jo Jo 

where / > was exploited and the asymptotic formula ^{r) = Cr'+^(1 + o(l)), r coming 
from equation ([3]) was used to estimate the integral. 

Since ipi{x) G L^(0,oo) and Q{x) G L^(0, oo) we infer that the m-function of the Sturm- 
Liouville operator ^[(^(x)] satisfies 

^ / (/ + l/2)^ \ ^ ^ ^ ka JUi/2ika)-t^n6iY{^,^,{ka) 

V J V/(0) c Ji+i/2{ka) -taii6iYi+i/2{ka)' 

From this we can build the potential through the formula relating the m-function to the spectral 
function and the constructive method discussed previously. 
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3.3 Deriving a moment problem 

With reference to the defining formula ([32D for F{x) we define a truncated version F{x) 



F{x) 



cos 



{VXx)da{X). 



(41) 



It turns out that the reconstruction of the F{x) function from the given m-function values is 
an inverse moment problem. Consider 



F{x)e^^^^)-dx 



^0 



C 



da{X)- 



' + +A 
(42) 

This integration can be performed with ease in terms of the m-function and by considering that 
on (— oo, 0) dp^{\) is concentrated at the bound states Ai, . . . , supported by S[Q{x)\. 



1 



- d(p^(A)-p°-°(A)) 
1 



+ 



dp^(A) 



1 



(43) 
(44) 



Here 



6, =p'^(A, + 0)-p'^(A,-0) (45) 

and mo(-) denotes the m-function associated with Q{x) = 0. The values of m(-) appearing in 
(1441) are determined by (HOj) . For a more detailed derivation consult [4]. 
Now we have the following problem for the truncated F{x): 



F{x)ey+-^)Ux = fi{c, h, 5u {A,}, {6,}), / = 0, 1, 
with the moments (assuming B bound states) 



(46) 



lj,{c,h,6i,{\i],{hi}) = / + 



J'i+i/2i^(^) - tan 6iY{^y^{ka) 



-1 



Ji+i/2{ka) - tan 6iYi+i/2{ka 



(47) 



which is essentially a moment problem. Note that unless there are no bound states the moments 
depend on the undetermined quantities {Aj} and {bi} associated to the bound state positions 
and norms. 



4 Solution method 

Depending on the number of bound states present in the auxiliary problem, different solution 
strategies are called for. As a consequence it is important to know or assess the number of 
bound states before solving the inverse problem which is discussed in the next section. 
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4.1 No bound states 

In this case we have a proper moment problem for F{x) = F{x): 



poo 

Jo 



(48) 



where the moments now do not depend on unknown quantities. 

To solve this moment problem we use the following expansion for F{x) 



N 



F{x) = y^c^e" 



(49) 



n=0 



Upon substitution into equation (HSj) and using + 1 fixed energy phase shifts as input data 
we get the following system of linear equations for the coefficients: 



N 



-cn + l + ^ 

n=0 ^ 



/ = 0, 1, . . . , iV. 



(50) 



Solving this system of linear equations yield the coefficients required to build F{x) and 
from that one can calculate the fixed energy potential essentially by solving the GL integral 
equation (130|) . Note that solving the system of equations is not a well-conditioned task as we 
must deal with a Hilbert-type matrix which is infamously badly conditioned. Therefore, from 
the numerical point of view it is of considerable value to see, that this matrix can be inverted 



explicitly. Our matrix, i.e. 



-cn+l-\ 



In 



is in fact a Cauchy matrix. The elements of the inverse 



of a general Cauchy matrix with elements Oj 



[X,, 



+ Vj) ^ are given by [H] 



Xm ~l~ Vi 
Xm Xj 



(51) 



In our case this implies 



N 



1 /, 1\ T-r I - cn' + ^ ^ I' - cn + I 



2 y -'- -'- cn — cn' /' — / 



n = 0,l,...,N. 



(52) 



Improvements of the solution method To further improve the solution method outlined 
above one can exploit two properties of the F{x) function, namely that 



F(0) = -h, 
lim F{x) = 0. 

a;— >oo 



(53) 
(54) 



The latter is a simple consequence of the Riemann-Lebesgue lemma, however the former needs 
more explanation. Let us write up 



da{\) = lim / da{\) = lim 

-oo J-oo 



/(A) _ p/^(_oo) 



TT 



lim 

A— s>oo 



-VA + p\-oo) -h + o(l) - p^-oo) - -VA 

TT TT 



-h (55) 

where equation fl2^ was used and continuity in A was supposed. This proves F{0) = —h. 
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Because of F(oo) = one can take 

Co = 0; (56) 

however, we note that this is not always the best course of action in practical scenarios (see 
Examples) . 

To incorporate the information -F(O) = —h there are several ways to choose from. For 
instance, one can prescribe the condition 

N 

5^c„ = -/i (57) 

n=0 

for the coefficients (which complicates the solution process: the matrix to be inverted is no 
longer of the Cauchy type). On the other hand, this will be very useful in the one bound state 
case (see later) where it will permit to determine the nonlinear parameter A in a linear way. 

Hausdorff moment problem Our problem can be viewed as a Hausdorff moment problem 
since with z = e'^ and ^{z) = —cz'^^^F{c\ogz) fHHl) takes the form 



/ z'd{z)dz = fii, / = 0,1,2,.... 
Jo 



(58) 



The Hausdorff moment problem is studied in the literature in detail concerning both mathe- 
matical properties and solution methods. For instance, an interesting stability result can be 
found in ^lOj whose corollary is the following theorem. It establishes an accuracy estimate of 
the inverse moment problem which is procedure independent. 

Theorem 1. Suppose the smoothness condition 

[\d'{z)-d'Niz)\'dz<E' <oo. (59) 

Then if the first + 1 moments of ^{z) and ^n{z) conincide, i.e. 

[ z^d{z)dz= [ z''dN{z)dz, k = 0,l,...N, (60) 
Jo Jo 

we have 

l'\^(')-^^(')\"d^<^(j^2- (61) 

Using Theorem [T] one can conclude that if the moments /i; are free of error, the difference 
between the approximated and the true F{x) functions in norm tends to as the number 
of moments is increased: 

|F(x) - FMl'dx < (62) 
with some C constant depending on the smoothness of F{x) — Fn{x), 

\F'{x)-F'j^{x)\^e-^dx<C, (63) 







and F]\f{x) is the approximation of the true F[x) using the first + 1 moments. 

It is interesting to consider the case when the input data is noisy. Using Theorem 1 of [TU] 
one can state the following: 

[ \F(x) - FN(x)\''dx < min | ^e^'^^^+i) + ^ : n = 0, 1, . . . A^l , (64) 
Jo " l|c| 4:{n + iy J 

where e"^ is the absolute square sum of the differences between the true and noisy moments. 
This result implies in particular that even if the number of phase shifts grows to infinity the 
recovery will not be complete when the data remains erroneous. 
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4.2 One bound state 

Supposing din]) we get the system of equations 



N 



— c 



-cn + / + 7j 

n=0 ^ 



+ 7 



/ = 0,1, 



iV + 1, 



(65) 



where fii denotes the Ith moment without the bound state contributions and A < and b > 
are the bound state parameters (the subscript 1 is omitted). Using the expansion fH9|) for F{x) 
we have 



N 



F{x) = bcosh.{V—Xx) + c^e" 



(66) 



n=0 



To get an exphcitly solvable system of equations treating a/— A as a parameter we subtract the 
term |e~^'^^ from the expansion for F{x), that is 



Fix) 



N 

E 

r),=0 



N 



-Xx 



Fix) = -e^ 



-Xx 



+E 



Cn.C 



(67) 



n=0 



and obtain 



-c6 



AT 



— c 



2((/ + i)+cv/^) _o 
through the elementary identity 

a 1 1 



■cn + I + 



1 1 



The explicit solution of the system of equations is given by 

. N+l 



1 = 0,1, 



.N + l 



(68) 



(69) 



1 



[I - cn + ^ + c6n-r 



^ n 



I — cn' + ^ + c6n'-iV—^ 



—cn' + c6n' -1 v—A + cn — c6n -i 



-A 



X 



n ''""'^!,^f"-'^ . (70) 



I and 5a,b is the Kronecker-delta. To determine a/— A 



where n = — 1 is also allowed, c_i 
we use a result concerning Cauchy matrices, namely that the sum of the coefficients (including 
6/2) is linear in \/—\ (see e.g. jl]. Lemma 5): 



N 



S = Cn = a a/— A + /3. 



(71) 



n=-l 



Using this fact the linear solution procedure is performed as follows. 



1. Calculate S for two arbitrarily chosen A values (Aqi and A02) solving f l68|) . Denoting the 
two sums 5*1 and S2, using the relation S = -F(O) = —h get V— A by 



-A 



iS2 + h)V^i-iSi + h)^ 



02 



(72) 



5*2 — 5*1 

2. Calculate the c„ coefficients from (168|) with the appropriate A parameter obtained in (i). 

3. Calculate g(r) through F(x) and the GL integral equation. 
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4.3 Multiple bound states 

For two or more bound states we propose a non-linear problem as follows. 

1. Solve the nonlinear system of equations 

B ^ N 

gstTTFr^T^ + S'-^^^^"" '^^ '''' 

B N 

5^&. + $^c„ = -/i (74) 

1=1 i=0 

with the index set L = [0, 1, . . . , + 2B — 1] in the variables Ai, A2, . . . , As, 61, 62, • • • , 

Co, Ci, . . . , Cn- 



B , N 

h 



2. Calculate F{x) by 

i=l n=0 

3. Solve the GL integral equation to obtain the potential. 

Numerically, it is worthwhile to initialize the variables for the nonlinear solver with the ones 
obtained supposing a vanishing potential q{r) = 0. This is because, roughly speaking, g(r) is 
expected to influence only the tail of the Q{x) potential and the bound states shall be close to 
those corresponding to g(r) = (see next section). 



5 Assessment of bound states 

We shall make observations based on the transformation formula for the potential f lT^ : 

Qix) = ^\qiae-^)-e], g(a) = 0. (76) 

No bound state is expected when q{r) > fc^ on < r < a, that is when Q{x) is positive 
everywhere, at least for h = 0. This is the case, for example, when we have a constant fixed- 
energy q{r) potential in the form 

g(r) = |^ Ok^. (77) 

10 for r > a. 

However, in practical applications it is a far more natural assumption that q{r) = already 
on b < r < a with some 6 > 0. In this case we have 

Q[x) = - (^^^ e"^H, < a; < -|c| log , (78) 

and only the tail of Q{x) is influenced by the fixed-energy potential g(r), e.g. Q{0) is solely 
determined by ka. For this reason with given k and a parameters by taking q{r) = we can try 
to calculate the approximate bound states or at least estimate their number. In what follows 
we show exact results for the constant g-potential, which case contains the zero g-potential case 
as well. 
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For Q{x) = — se~^*^ with s,t > the Sturm Liouville equation can be solved exphcitly. 
This is an easy exercise and will not be detailed here, only the result is given. The differential 
equation 

- ^"(x) - se-2*Xx) = \ip{x) (79) 

is solved by 



t 



+ C2 J_, 



t 



'tx 



(80) 



where Bessel functions generally of complex orders have appeared. For ImA/A > the l?'- 
solution (if t > 0) is 



-tx 



Since 



J 



±is/\lt 1 ^ 



— e = h{x)x 



=Flmv^/i _|_ ^j^^=Flmv^/i^ 



X 0, 



(81) 
(82) 



where h[x) is a bounded function. Then the m-function is 



m(A) 



J'-.^.,t[^ 



J 



(83) 



Using the Stieltjes inversion, equation dp^{X) can be found. For A > trivially 



rfp'^(A) = -Im 

71 



J-iVx/tiVs/t) 



dX, 



A > 0. 



(84) 



For A < the measure is concentrated to points which are the bound states. They are 
located at the eigenvalues of the operator where ip G L^(0, 00). Starting from ip{0) = 1 and 
ip'{0) = h one can calculate the coefficient of the diverging solution to be 



TTWS 



2t 



(85) 



where some elementary properties of the Bessel functions were used [8]. Ci needs to be zero 
thus the bound states are located at A's that satisfy 



(86) 



Note that the number of bound states is greater than zero but finite. 

Using the theory of residues the height of the step in p(A) at the bound states can be 
obtained to be 



p(Ao + 0) - p(Ao - 0) = 2tv^ 



J, 



-\o/t 



(87) 



where the superscript (n, m) means derivation with respect to order n times and derivation 
with respect to the variable m times. 

In case of /i = we have a particularly simple scenario at hand. For definiteness let 



n^ = k^- g(0) 



(88) 
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Table 1: The first few sectors of definite bound state numbers for a constant q{r) potential at 
h = 0. 



0< Ka< 3.8317 
3.8318 <Ka< 7.0156 
7.0156 < Ka< 10.174 
10.174 <Ka< 13.324 



one bound state 
two bound states 
three bound states 
four bound states 
etc. 



(g(0) being the value of the constant potential at the origin, and suppose for simplicity that 
g(0) < fc^). The potential Q{x) = — ^^^^e~Tt with h = has bound states Aj at 

J'\c\V^S^^) = 0- (89) 

From [8] we infer, denoting the nth root of J^(a;) by „, that „ < j'^+e,n n = 1,2, . . . holds 
with J'qi = 0. Moreover, is continuous in /i [8J; therefore, we find that for na > there 
always exits at least one bound state. From the above fact it also follows that the number 
of zeros of Jq{x) on < x < equals the number of bound states of Q{x) = — -^^^-e"!^. 
Alternatively, as Jq{x) = — Ji(x), the number of bound states increase at the zeros of Ji{x), 
i.e. at ji^n, n = 1,2, . . .. Some numerical results are listed in Table [T] For a general potential 
g(r) the data listed in Table [T] can still be relevant if the potential is shallow comparad to k"^ 
and only the tail of Q{x) is influenced by q{r). Note that the value of the parameter c does not 
affect the number of bound states (only their positions) as it is only a scale parameter of the 
function of \/— A whose zeros give the bound states. 

Allowing h ^ gives rise to the more involved condition fl86|) for the bound states. It is 
possible then for some h ^ that one has one bound state while for h = two of them. This 
has the favorable consequence of reducing a nonlinear problem to a linear one. Without giving 
an exhaustive treatment of the situation we show the following illustrative result for bound 
state reduction. 

Lemma 1. For 3.83 ~ < 'f^' < jo,2 ~ 5.52 the two bound states present at h = can be 
reduced to one by varying h. 

Proof. In figure [1] the two Bessel-type functions entering condition fl86|l are depicted. We will 
show that by varying h (on the figure through R), disregarding an overall sign, one can always 
get the same kind of graphs as the ones on the figure. 

By using the fact that the zeros of Jo{x) and Jq{x) interlace we infer Jo^na) ^ and 
J^lna) 7^ 0. The only thing that remains to be shown is that the distribution of the zeros are 
as depicted. Since < na < jo,2 for J'f,{na) only the first two and for Jk{na) only the first 
zeros can enter our considerations. Then the interlacing relation j'^ ^ < ja,i < j'a 2 translates to 
ki < k2 < ks if I = na, jk2,i = i^cl and j^^ 2 = This completes the proof of the lemma. □ 

6 Examples 

6.1 Reconstruction of constant potentials 

First we reconstruct a potential 

fl.2 for r < 2 = a 

. , (90) 
10 for r > 2 = a 
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Figure 1: Julna)' (full line) and RJ^{K.a) (dashed line) as functions of the order {R is a real 
number satisfying sgn(jQ(Ka)) = sgn(i?Jo('f2)) and |i?Jo(Ka)| > |jQ(Ka)|). While at /i = the 
zeros of the derivative functions are the bound states, at /i 7^ the intersection of the two 
graphs give them. 

at = 1 that generates no bound states in the auxiliary problem (because Q[x) > 0). 

Figure!^ shows the quality of the reconstructions when 5, 10, 20 and 40 input phase shifts 
are used with the standard choice of parameters (c = —1, h = 0). As expected we see that the 
quality of the inversion procedure gets better as the number of input data is growing. 

We may let the parameters c and h change from the standard values and prescribe the 
smoothness of the inverse potential q{r). This can be done by monitoring the quantity 



at each values of the parameters c and h where the lower limit Tq is introduced in order to 
exclude (possible pathological) singularity at the origin. Note that we used rg = 0.05 in all 
the examples shown. By requiring s to be small is equivalent to exclude nonphysical potentials 
(see also [H]). On the other hand, this requirement also means that we encounter smooth 
F-functions too, in agreement with Theorem [TJ 

Figure [3] shows that choosing a modest number = 10 of input data we also get excellent 
reconstruction when s is small. The optimal choice of parameters has proved to be c = —0.3 
and h = -0.15 (Hd)). 

Several other techniques can also be employed to find optimal results. One is to retain the 
constant term cq in the expansion fH9|) . The value of cq (i.e. the departure from zero) serves 
also as an overall indicator of the accuracy of the procedure. Another trick is to reconstruct 
the potential using the one bound state procedure, that is by retaining the c_i = 6/2 term and 
the associated nonlinear parameter a/— A (a possible spurious bound state). This procedure 
also affects beneficently the numerics and can be applied both in the zero and one bound state 
cases. 

In Table [2] we list the input phase shifts 61 and the intermediate data (moments fii, coeffi- 
cients Cn and the value \/— A) of the inverse calculation with c = —0.3 and h = —0.5 (where 
s = 0.0004). As we see the procedure yields a/— A = —1.4447 for the "bound state" parame- 
ter, a negative value which clearly indicates that there is no true bound state in the auxiliary 
problem. 

^For brevity of the figure captions we introduce tlie function Ha{x) being a step function: Ha{x) = 1 for 
X < a and Ha{x) — for x > a. 
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Figure 2: Reconstructions with different number of input phase shifts of the constant potential 
q{r) = 1.2 H2ir). (c = -1, h = 0.) 



Now we proceed to reconstruct the potential 

f 0.8 for r < a 
qir) = \ - (92) 

I U for r > a 

at = 1 that generates one auxiliary bound state with a = 2 and two auxiliary bound states 
with a = 11 (choosing h = 0). 

Figure im^a)-(b) shows that the inversion procedure yields a better potential for the smaller s 
value also in this one (auxiliary) bound state case. The deviation of the inverse potential from 
the original one in the region 0.5 < r < 2 is of the order of 0.00010 for s = 0.049 (at c = — 1, 
h = 0) and 0.00005 for s = 0.0022 (at c = -0.5, h = -0.65). 

Figure 111(c)- (d) shows two reconstructions in the case when two bound states exist in the 
auxiliary problem. In the first case (figure Hl^c)) the inverse potential is obtained by the use 
of the two bound state method with parameters c = —1.5, h = and finding the nonlinear 
(bound state) parameters at Ai ~ —0.48, A2 ~ —2.43 (starting from the trial values of Ai ~ 
— 0.8, A2 ~ —2.3). In the second case (figure |l](d)), the inverse potential has been calculated 
by using the one bound state procedure with the parameters choice c = —1, h = 5 (providing 
for the bound state parameter the value A = —2.41). We see that both methods yield similar 
results. However, while the two bound state calculation needs an a priori guess about the 
positions of the auxiliary bound states, the one bound state method does not suffer from such 
an ambiguity and the calculation can be performed in the linear regime. One can check that 
in this case na ^ 4.92, thus we are in the domain where the number of bound states can be 
reduced (Lemma 1) by varying h. 
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(a) c = —1, h — (standard choice). s=12. 




(b) c= -l,h^ -0.15. s=:15. 




1.00 




(C)C: 



-0.30, h = 0. s=0.15. 
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(d) c = -0.30, h = -0.15. s=0.011. 



Figure 3: Reconstructions with different choices of the parameters of the constant potential 
q{r) = 1.2H2{r) from 11 phase shifts with no bound states in the auxihary problem. 



6.2 Reconstruction of potentials with different shapes 

Figure [5t^a)-(b) shows the reconstruction of the two potentials given by 



gGauss(?") = -4e qwsir) = - 



4 

~ r-0.5 

1 + e 0.1 



(93) 



The reconstructions have been carried out at k = 1.5 with 7 input phase shifts given with a 
precision of 4 digits for the Gauss potential O'caussl'"), and with 4 input phase shifts given with 
a precision of only 2 digits for the Woods-Saxon shaped potential gws('^)- In both cases the one 
bound state approximation procedure is applicable. For the Gauss potential the reconstruction 
with the parameters a = 1.5, c = —0.74, and h = resulted in the bound state position at 
A = —3.36. Note that this figure agrees with the exact bound state position (when calculated 
from the known auxiliary potential Q{x)) and only slightly differs form the assessment value of 
A = —3.22 (calculated from equation ( 15^ ). For the WS potential the reconstruction was carried 
out with the parameters a = 2,c = — 1.25,/i = and resulted in the bound state position at 
A = —2.46 to be compared with the assessed value of A = —2.44. 



6.3 Inverse potentials from experimental phase shifts 
6.3.1 e- Ar atom scattering at E = 12 eV (0.4412 au). 

The experimental phase shifts derived by Williams [12j from e — Ar scattering experiment at 
Ec.m. = 12 eV and the inversion parameters a and c are listed in the first line of Table |3l At 
c = —3.7 and h = 1.9 the bound state parameter was found to be A = —0.44. The resulted 
potential of the HA calculation is shown by the continuous line in Fig. El^a) and compared with 
that obtained by the modified Newton-Sabatier (mNS) inversion method [13] (dashed curve). 
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Table 2: Input phases 6i and the intermediate quantities of the inversion procedure of the 
potential q{r) = 1.2 H2{r): moments Hi, coefficients Cn and the nonlinear parameter (i.e. 'bound 
state'). 



/ 


Si 










n 








— 1 




1 

X 


—0 99fi4 




—0 0043 


n 


_i_n noo9 


2 


-0.0471 




0.0151 


1 


-0.0954 


3 


-0.0037 




0.0180 


2 


+7.2085 


4 


-0.0001 




0.0176 


3 


-0.1510 


5 


-5.0 X 10- 


-6 


0.0164 


4 


+0.2657 


6 


-1.1 X 10- 


-7 


0.0151 


5 


-0.2968 


7 


-1.8 X 10- 


-9 


0.0139 


6 


+0.2027 


8 


-2.2 X 10- 


11 


0.0129 


7 


-0.2022 


9 


-2.3 X 10- 


13 


0.0119 


8 


+0.0258 


10 


-1.9 X 10" 


■15 


0.0111 


9 


+0.0092 



^ Nonlinear "bound state" parameter: \/—\ = —1.4447. 

Table 3: Input phase shifts Si derived from e — Ar and n — a scattering experiments performed 
at one and three centre of mass energies, respectively. Corresponding inversion parameters c 
and a. h = 1.9 was used in case of the e — Ar inversion while h = was taken for the three 
n — a inversions. 





5o 


5i 


^2 ^3 


a 


c 


12 (eV) 


-1.218 


-0.626 


1.191 0.118 


3.9 


-3.70 


9.6 (MeV) 


1.763 


1.553 


0.028 


3.9 


-4.25 


12.8 (MeV) 


1.676 


1.466 


0.066 


3.4 


-2.96 


16.0 (MeV) 


1.588 


1.396 


0.117 


3.3 


-2.37 



Both potentials compare well showing the usability of the present method. The e — Ar potential 
has an attractive part with a minimum value of about —2.8 au at a distance of r ~ 1.2 au. At 
smaller distances the potential is of repulsive nature which can be interpreted as a manifestation 
of the Pauli-principle. Further details can be found in Ref. |14j . 

6.3.2 n — a particle scattering. 

Input data at energies Ec.m. = 9.6, 12.8, and 16.0 MeV are listed in the lower three lines of 
Table [3] together with the inversion parameters a and c. The input phase shifts are taken from 
the comprehensive analysis of the n — a scattering presented by Ali et al [15] . Because of the 
spin-orbit coupling both spin-up 6^ and spin-down phase shifts contribute to the scattering 
amplitude at each partial wave. In case of weak spin-orbit coupling (which is assumed) the 
combined phase shifts 

Si = ^{(i + + (94) 

are characteristic of the underlying central potential [16], and these data are used as input for 
the HA procedure. 

The resulting three HA potentials are exhibited in Fig. El^b). As we see they offer a similar 
physical interpretation for the scattering process as before in the case of electron-inert gas 
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(a) c = -1.0, /i = 0, A = -0.105, s = 0.049. 



0.810 




0.790 



(c) c = -1.5, /i = 0, Ai = -0.48, A2 = -2.43, 
s = 1.28, two bound states. 
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(b) c = -0.5, h = -0.65, A = -1.39, s 
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(d) c = -1.5, /i = 5, A = -2.41, s = 0.11, one 
bound state. 



Figure 4: (a)-(b) Reconstructions of the constant potential g(r) = 0.8H2{r) from 11 phase 
shifts and one bound state in the auxihary problem at h = with different choices of the 
parameters, (c)-(d) Reconstructions of the potential g(r) = 0.8ifii(r) using (c) the two bound 
state formulation {h = 0) and (d) the one bound state method {h = 5). The constant term was 
present in ( H9|) during the calculation. 



atom collision: the approaching colliding partners, neutron and alpha particle are attracting 
each other when entering the domain of nuclear forces. The attraction is culminating around 
the a— particle surface between r f» 1.1 — 1.2 fm, reaching a strength of potential energy between 
—50 and —52 MeV. When the colliding partners are merging the Pauli repulsion (originating 
from the fermionic exchange) takes into effect and this is overcome by the nucleonic soft core 
repulsion at very small distances at r ^ 0.1 — 0.2 fm. Because the resulting inversion potentials 
are also very similar at these different energies between 9 — 16 MeV, we may have found the 
energy-independent potential responsible for the scattering data which is always the desirable 
goal of any fixed energy inversion procedure. 



7 Conclusions 

We have surveyed, developed further and applied the constructive inverse scattering method of 
Horvath and Apagyi (HA) The inverse problem consists of finding a scattering potential 
of finite support in the radial Schrodinger equation from a finite number of scattering phase 
shifts given at one fixed energy. The solution procedure consists of solving an auxiliary inverse 
spectral problem of the classical Sturm-Liouville equation whose spectral data is determined by 
the phase shifts of the fixed energy inverse scattering problem. The auxiliary inverse spectral 
problem is equivalent to a moment problem. 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 0.5 1.0 1.5 2.0 

r r 

(a) Gauss potential, s = 4.9. (b) Woods-Saxon potential, s — 7.4. 

Figure 5: Reconstruction of potentials with different shapes, (a) Gaussian shape from 7 phase 
shifts, (b) Woods-Saxon shape from 4 phase shifts. The original potential is depicted with 
dashed line while the reconstruction is illustrated with solid line. For parameters and accuracy 
see the text. 




Figure 6: Construction of potentials from experimental phase shifts listed in Table |3l (a) HA 
(continuous line) and mNS (dashed line) inverse potentials for e — Ar scattering at 12 eV. (b) 
Inverse potentials for n — a scattering at three cm. energies of 9.6, 12.8, and 16.0 MeV. 



The HA method has been developed further in that two hidden parameters of the theory have 
been disclosed and enabled to vary. One of these parameters is a scale parameter c appearing 
in the Liouville transformation. It plays an important role in the solution of the moment 
problem. The other parameter is an initial value parameter h which enters the Gel'fand-Levitan 
constructive inversion scheme. The hidden (also called standard) values of these parameters 
have been c = — 1 and h = 0. By making these parameters free the HA method proves to be 
applicable to cases where only a limited precision of input data is expected as in the case of 
the inversion of phase shifts derived from measured cross sections. 

Depending on the number of bound states present in the auxiliary inverse spectral problem, 
various solution methods of the moment problem have been presented and applied to calculate 
the potential from a set of given phase shifts which can be either calculated theoretically if 
the underlying potential is known (reconstruction procedure) or derived, e.g. from collision 
experiments if the potential is not known (construction procedure). Of course, an inverse 
method aims at solving the latter task, and the prior one may serve for testing or developing 
further the procedure. 

The examples illustrate that a dramatic improvement of the HA inversion method can 
be achieved by the proper adjustment of the parameters c and h. We may, for example. 
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reduce the number of bound states of the auxihary inverse spectral problem by one. Because 
the most sensitive part of the method is the possible presence of auxiliary bound states we 
have established a procedure to assess their number. This procedure is based on the free 
motion because only the asymptotical part of the auxiliary potential is influenced by the wanted 
fixed energy potential. Alternatively, one may also use other (less sophisticated) fixed energy 
inverse scattering methods, e.g. the modified Newton-Sabatier procedure [T3|, to assess the 
number of auxiliary bound states. By using an independent phase equivalent method, one can 
simultaneously check the inverse potential provided by the HA procedure. 

Finding an optimal choice of the parameters a, c, and h can proceed through satisfying phys- 
ical arguments. For example, one can obtain a cheap calculation (that is using small number 
of input data) via prescribing a maximal smoothness for the inverse potential. With such a 
prescription we have reproduced known potentials. Another (or simultaneous) prescription can 
be to check the inverse potential whether it fulfills the approximate relation V{a) ~ 0. With 
this prescription we have got the inverse potentials equivalent to the measured phase shifts of 
e — Ar and n — a scattering experiments. Note that by making the parameters c and h flexible 
one may also avoid the appearance of non-physical (e.g. singular) inverse potentials [TT] which 
can arise e.g. when the solution of the GL equation is not unique. 

The HA method can be extended into various directions. A natural extension is to develop 
the theory to handle complex phase shifts which describe inelastic processes. Another exten- 
sion would be to efficiently treat Coulombic processes, i.e. charged particle scattering. An 
interesting and important extension would be an alternative formulation where the auxiliary 
spectral problem is solved by the Marchenko integral equation [T7]. The latter development is 
in progress. 
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